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ABSTRACT 

Numerous cosmological simulations have been performed to study the formation of the 
first objects. We present the results of high resolution 3-D cosmological simulations 
of primordial objects formation using the adaptive mesh refinement code FLASH by 
including in an approximate manner the radiative transfer effects of Lyman alpha 
photons. We compare the results of a Lyman alpha trapping case inside gas clouds 
with atomic and molecular hydrogen cooling cases. The principal objective of this 
research is to follow the collapse of a zero metallicity halo with an effective equation 
of state (that accounts for the trapping) and to explore the fate of a halo in each of the 
three cases, specifically, the impact of thermodynamics on fragmentation of halos. Our 
results show that in the case of Lyman alpha trapping, fragmentation is halted and a 
massive object is formed at the center of a halo. The temperature of the gas remains 
well above 10 4 K and the halo is not able to fragment to stellar masses. In the atomic 
cooling case, gas collapses into one or two massive clumps in contrast to the Lyman 
alpha trapping case. For the molecular hydrogen cooling case, gas cools efficiently and 
fragments. The formation of massive primordial objects is thus strongly dependent 
on the thermodynamics of the gas. A salient feature of our results is that for the 
formation of massive objects, e.g. intermediate mass black holes, feedback effects are 
not required to suppress H2 cooling, as molecular hydrogen is collisionally dissociated 
at temperatures higher than 10 4 K as a consequence of Lyman alpha trapping. 
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1 INTRODUCTION 

The formation of primordial objects is not yet very well un- 
derstood. A lot of progress has been made in analytical and 
numerical approaches but still further investigation is re- 
quired. The formation of the first objects depends on the 
fragmentation process and its relation to ther modynamical 
properties of the gas ([Ciardi & Ferrarg] I2005T ). In order to 
understand their formation, it is crucial to investigate the 
impact of thermodynamics on primordial halos which define 
the mass scale of first objects. It is generally believed that 
first stars are formed in minihalos of ~ 1O 6 M0 at redshift 
z rsj 20 — 30 while qu asars are formed in massive halos at 
z - 10 (|Brommll200l ). 

Numerical simulations performed to 
mation and fragmentation of 



study the for- 
clouds 
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in ~ 10 M halos if gas is condensed by H2 cooling. It 



has been seen in these simulations that first stars are 
massive due to inefficient cooling (in the absence of metals 
and dust) a nd ar e formed in isolation. Recent work by 
IClark et al.1 (|2010h shows that primordial gas is susceptible 
to fragmentation even with small amount of turbulence. 
For metal free gas, the only efficient coolants that can bring 
down the temperature below ~ 10 4 K are H2 a nd HP. Gas 
will n ot fragment to form stars in their absence (|Abel et al.l 
l2000h . Feedback processe s can suppress the formation of 
molecular hydrogen ([Haiman et al.ll200Cft . 

Cooling due to molecular hydrogen in halos having 
virial temperatures > 10 4 K can be suppressed through 
Lyman- Werner radiation emitted by a nearby star forming 
galaxy. Such hal os therefore could be suitable for the growth 
of black holes ([Diikstra et al.l [2008) . Early fragmentation 
will not occur when H2 cooling is suppressed either due to 
external or inter nal sources and will favor the formation of 
massive objects (lOh fe Haimanll2002l : iBromm fc Loebll2003l : 
iBegelman et al. 1 120061 ). In the absence of molecular hydro- 
gen cooling, gas can only cool through atomic line radia- 
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tion. I Wise et alJ |2008 ) performed cosmological simulations 
to study the collapse of protogalactic gas clouds cooled by 
only atomic lines and found that an object of 10 5 M© is 
formed in the center of a metal free halo. They also deter- 
mined that angular momentum tra nsfer played a vital role 
in the formation of a massive object. iBromm fc Loebl <l20Q3h 
used SPH simulations to study the evolution of metal free 
gas in the absence of H2 cooling. Their results show that col- 
lapse is isothermal and fragmentation is not very efficient. 
As a result one or two ma ssive clumps of 10 6 Mq are formed 
containing > 10% of gas. iRegan & Haehneltl (2009b : a) car- 
ried out cosmological simulations and followed the dynamic 
evolution of metal free halos cooled by atomic line radiation 
and found that isothermal collapse leads to the formation of 
a self gravitating stable disk at the center of a halo. 

Massive black ho les can only form in atomic cooling ha- 
los ([Abel et al.ll2000h with virial temperatures greater than 



10 K. Self- gravitating gas can form black holes through 
direct collapse if it is cooled down to a small fract ion of 
the virial temperature dBegelman fe Shlosm an 2009). Var- 
ious other w ays have been su ggested for the formation of 
black holes (iHaimanl l2004bl l 3; IPortegies Zwart et al. _[2004 ; 
Devecchi Volonteril I2OO9I : iReesI Il978l : iDiorgovski et all 



2008). Gas in primordial halos can cool through atomic 
line radiation and collapses isothermally to form a 10 4 — 
10 6 Mp object if the te mperature remains as high as 10 4 
K (|Spaans &; Silkl l2006h . These authors found that trap- 
ping Lyman alpha photons stiffens the equation of state and 
raises the temperature above 10 4 K. The fraction of molecu- 
lar hydrogen is almost zero at such higher temperatures. 

Trapping Lyman alpha photons increases the Jeans 
mass and fragmentation is likely to be inhibited. The im- 
pact of Lyman a lpha pressure on accre tion of first stars has 
been studied bv iMcKee fc Tanl (|2QQ8h . They found that it 
can not effect overall infall except in polar regions of 20-30 
M0 stars. The key assumption in most models for the for- 
mation of massive objects through direct collapse is that gas 
must avoid fragmentation otherwise it might form low mass 
Pop III stars. Hence, que nching the formation of H2 will be 
necessary. The model by S paans fe ZsH3 <l2006h for the for- 
mation of pregalactic black holes requires no explicit mecha- 
nisms for the destruction of H2 as trapping of Lyman alpha 
photo ns keeps the gas H2 free, but also see ISchleicher et al.l 
(l2010h . Motivated by their model, we have performed 3-D 
cosmological simulations that include trapping effects of Ly- 
man alpha photons and compared our results with previous 
work (atomic and H2 cooling cases). The only assumption 
for this model is that the halo is metal free. The metal dis- 
tribution is generally inhomogeneous and metal enrichment 
has not proceeded very far at higher redshift s. Pockets of 
meta l free gas can exist thus until redshift six ([Trent i et al.l 
l2009h . 

The primary aim of this research is to follow the col- 
lapse of a (close-to) zero metallicity halo and to explore, 
whether it fragments into multiple system or forms a mas- 
sive object with out fragmenting. We carry out high resolu- 
tion cosmological simulations that have a dynamic range of 
~ 2x 10 5 in linear scale and study the impact of thermody- 
namics on the collapse of primordial gas. We present in this 
paper three sets of high resolution cosmological simulations 
with different physics ingredients. In the first case, we add 
atomic cooling assuming that there is a UV background that 



suppresses molecular hydrogen formation. In the second set 
of simulations we include cooling due to molecular hydro- 
gen to see how it affects the process of collapse. In the third 
suite of simulations, we assess the influence of Lyman al- 
pha photon trapping on the formation of primordial objects. 
We compare the results of the Lyman alpha trapping case 
with atomic and molecular hydrogen cooling cases. For the 
case of Lyman alpha trapping, we assume that no molecular 
hydrogen formation takes place at temperatures above the 
collision break-up energy of H2. Also, even small amounts 
of dust, if present, can absorb and re-e mit the Lyman al- 
pha p hotons at much lower frequencies ([Haiman fe Spaansl 
1999). So, we assume zero dust abundance. 

Our paper is organized as follows. In section 2, we de- 
scribe the simulation setup and numerical methods used in 
this paper. In section 3, we discuss different heating and 
cooling mechanisms that we included in our simulations. Fi- 
nally, in section 4, we discuss the results obtained with our 
simulations and present our conclusions. 



2 SIMULATIONS 



We are using the code FLASH3 ([Frvxell et al.ll2000h . FLASH 
is a module based, Eulerian AMR, parallel simulation code 
which can solve a broad range of astrophysics problems. 
The Message- Passing Interface (MPI) library is used to get 
portability and scalability on many different parallel sys- 
tems. It makes use of the PARAMESH library to handle a 
block-structured adaptive grid, to add higher resolution de- 
pending upon user's defined refinement criteria. For hydro- 
dynamic calculations we use the directionally split pi ecewise- 
parabolic method fPPM. IColella fc Woodward! 19841 ) and in- 
side shock waves we use the HLLE solver. PPM is favored 
above other methods due to its higher order accuracy and 
its ability to sharpen shocks and contact discontinuities as 
it does not use artificial viscosity. In order to calculate the 
evolution of dark matter, FLASH uses an N-body particle 
mesh technique and couples it to the hydrodynamics. The 
gravity solver used in FLAS H is the oct-tr ee based multigrid 
Poisson solver developed bv iRickerl (| 20081 ). 

We take a three-dimensional box of 10 Mpc (co-moving) 
on each side. In order to simulate the evolution of dark mat- 
ter, we are initializing 2.6 x 10 6 part icles. We use t he CO S- 
MICS (grafic) package developed bv iBertschingerl ([ 19951 ) to 
introduce Gaussian random field cosmological initial condi- 
tions. As motivated by inflation, the density perturbations 
are assumed to be isotropic, homogeneous and Gaussian. 
We use standard parameters from the WMAP 5-year data 
(Q m = 0.2581, H = 72kms" 1 Mpc"\ Q h = 0.0441) with 
a value of as = 0.8. We use periodic boundary conditions 
both for gravity and hydrodynamics. 

We run the COSMICS package to create Gaussian ran- 
dom initial conditions and select the set of initial conditions 
having the highest density peak. We find the location of the 
highest density peak and shift it to the center of the box. 
In the region of IMpc around the center of the box, we ini- 
tialize our simulation with a high level of resolution, which 
corresponds to 512 3 grid cells. The rest of box is set to a 
resolution of 128 3 cells. We start the simulation at redshifts 
ranging from z=50-70 for different sets of Gaussian random 
field initial conditions with higher and lower density peaks 
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halos. We select a massive halo and zoom in with the FLASH 
refinement criteria. As adding higher levels of uniform re- 
finement within the region of interest becomes a computa- 
tionally demanding task, we exploit the AMR technique to 
add dynamic levels of refinement. We add 8 additional dy- 
namic levels of refinement using specific density threshold 
criteria. This gives us 15 levels of refinement in total. The 
effective resolution in the region of interest is ~ 2 x 10 5 in 
linear scale (i.e. dynamic range). The virial mass of the sys- 
tem is ~ 3 x 1O 9 M0 and the virial temperature is ~ 10 5 K 
(T vir = 10 4 x (M v i r /M ) 3/2 (l + z vir /lO)). The collapse hap- 
pens at redshift 10. In the case of molecular hydrogen cool- 
ing, the Jeans mass reduces as gas cools to lower tempera- 
tures and collapes to higher densities compared to the other 
simulations, so we add an additional dynamic refinement 
level to resolve the Jeans mass (maximum refinement in this 
is 16 levels). Our refinement criteria are good enough to re- 
solve shock waves and they also ensure the fulfillment of 
the Truelove criteria stating that the Jeans length sh ould 
be resolved by at least four cells ([Truelove et al.|[l997h . In 
this way, we avoid artificial fragmentation caused by density 
perturbations created during discretization of the grid. The 
top level hierarchy of refinement levels helps in even greater 
precision. 



3 PHYSICS INCLUDED 
3.1 Cooling 

Cooling is an important process in the formation of struc- 
tures, like stars and galaxies. It is crucial to understand the 
physics of cooling in primordial gas in order to study the 
formation of first objects. Adequate amounts of gas has to 
gather in dark matter halos to form structures. Cooling be- 
comes efficient when the cooling time t CO oi = 3nkT/2A be- 
comes less than the free-fall time tg = y37r/32Gp, where 
n is the number density, p is the mass density and A is the 
cooling rate. Under this condition gas cools, collapses and 
can fragment. Mass scales of fragments are determined by 
the Jeans mass Mj w 1O 6 M (T/1O 4 K) 3/2 (for typical densi- 
ties of 10 4 cm -3 ) which depends on the thermodynamics of 
the gas. Gas cannot cool below 8000 K if it is cooled by only 
atomic line radiation. Consequently, the Jeans mass remains 
> 1O 5 M0 and fragmentation to lower mass scales is inhib- 
ited. In order to induce fragmentation of lower mass scales 
gas should be cooled down to lower temperatures. Introduc- 
tion of molecular hydrogen and HD cooling can lower the 
gas temperature down to <~ 100 K and reduces the Jeans 
mass to rsj 1O 3 M0. Fragmentation at lower than the Jeans 
scale is prohibited as gas becomes pressure supported. 

The net cooling function of a gas is a combination 
of processes, like line radiation, photo-ionization, compton 
scattering, bremsstrahlung, radiative recombination, etc. 
The cooling function strongly depends on the composition 
of the gas. The primordial composition of the gas is dic- 
tated by big bang nucleosynthesis. We take a primordial 
composition for the gas with 75 percent hydrogen and 25 
percent helium by mass. In metal- free gas the dominant 
cooling process is atomic line cooling. As primordial gas is 
metal- free, we are implementing a zero metallicity cooling 
function bv lSutherland &; Dopital (| 19931 ). At higher temper- 




Temperature [K] 

Figure 1. Cooling function for zero metallicity. The dashed line 
shows the zero metal cooling function while the dash-dot line 
depicts the contribution due to molecular hydrogen. 



atures of 10 6 — 10 8 K cooling is due to free- free transitions. 
For temperatures between 10 4 — 10 6 K cooling is mainly due 
to HI and Hell Lyman alpha. 

Evolution of primordial density fluctuations is con- 
trolled by the ability of gas to cool down to low temper- 
atures. The temperatu re of gas remains ~ 8000 K in the 
absence of H2 and HD ([Galli fe Pallalll998t ) due to high ex- 
citation energies of primordial gas atoms. In the absence of 
metals, H 2 cooling is crucial to form stars in T v i r > 10 4 K 
halos. Vibrational and rotational modes of H2 can be ex- 
cited at low temperatures <~ 500 — 1000K. Trace amounts 
of molecular hydrogen can form via gas phase rea ctions in 
the post-recombination era (Peeb les fc Dickelll968h . We use 
the lGalli fc Pallal (199S|) cooling function to include cooling 
due to molecular hydrogen. We take a universal abundance 
of H 2 , i.e., 10" 3 . The zero metal cooling function is shown 
in figure [T] details can be found in the above mentioned 
articles. 



3.2 Variation in the equation of state 

The equation of state plays a vital role in structure forma- 
tion. Fragmentation of gas clouds is str ongly depende nt on 
the polytropic equation of state (EOS) ([Li et al.ll2003h . The 
relation of the polytropic index of the EOS (7) to logarith- 
mic derivatives of heating and c ooling function impli citly in- 
cludes radiative transfer effects (jSoaans fc Silkl2QQ0h . Under 
conditions of thermal equilibrium, by balancing the heating 
and cooling terms, the polytropic exponent which measures 
the c ompressibility of gas can be written as ([Scalo &; Biswas! 
l2Q02h : 



7=1 + 



dlogT 



(i) 



dlogp 

ISpaans fc Silkl (|2006h have studied the effect of trapping Ly- 
man alpha photons on the polytropic EOS for primordial 
halos. They found that trapping Lyman alpha photons stiff- 
ens the polytropic exponent significantly above unity. Con- 
sequently, the Jeans mass is boosted and fragmentation is 
halted. 

As the collapse proceeds gas becomes optically thick. 
The escape probability of Lyman alpha photons, which per- 
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form a random walk both in space and frequency, is reduced 
and ultimately photons are trapped inside the cloud. If the 
photon travel time becomes larger than the free-fall time, 
the radiative energy is re-deposited inside the gas and con- 
sequently cooling is suppressed and the temperature of the 
gas becomes ^ 10 4 K. The formation of molecular hydrogen 
will not take place as it is col lisionally dissociated at tem- 
peratures higher than 8000 K. iRees have 
studied the influence of trapping radiation on collapsing gas 
clouds and concluded that trapping of Lyman alpha pho- 
tons can inhibit fragmentation. They also mentioned that 
trapped radiation would exert pressure which is half of the 
pressure required to support the cloud and will not halt the 
overall collapse. Radiation pressure produced by trapping 
of Lyman alpha photons c an considerably con strain the ef- 
ficiency of star formation (|Oh fc Haimanll2002h . 

At optical depths To > 10 7 the photon escape time t p h 
becomes larger than tff. Due to the weak dependence of t p h 
on the number density (t p h oc n -1 / 9 ), the photon escape 
time remains longer than the free fall ti me (t p h oc n _1//2 ) 
during the collapse (jSoaans Silkl l2006h . The number of 
photons that is trapped at line center thus increases and 
the emission of Lyman alpha photons from the core tends to 
zero. Consequently, trapping becomes more effective during 
the collapse. The scattering of Lyman alpha photons during 
the collapse can introduce extra pressure which slows down 
the process of collapse. This changes the density of trapped 
Lyman alpha photons and induces fluctuations in the EOS 
which are not well captured in our case. A temporary soften- 
ing of the EOS would favor fragmentation and our approach 
therefore gives an upper limit on the effect of Lyman alpha 
trapping. The minimum hydrogen column density required 
for the trapping is Nh = 10 21 cm -2 . 

By trapping Lyman alpha photons the polytropic index 
7 increases from 1.01 to 1.50 for hydr ogen number densities 
of 1 to 10 5 cm- 3 (|Spaans fe Silkl I2006T ). We found that high 
columns in our case sti ffen the EOS a t lowe r density val- 
ues than mentioned by ISpaans & Silkl {2006). Due to high 
columns, gas becomes opaque at lower number densities 
and causes more efficient trap ping of Lyman alpha photons. 
Therefore, we have modified ISpaans & Silkl (|2006h criteria 
for the change in polytropic index to include high column 
opacity effects. Their model is based on analytical results 
and does not take into account the resolution constraints. 
We use the formula given below to implement the change in 
the polytropic exponent of the EOS 7, 



+ T8 Bn l 



7/18 



7/18 ' 



(2) 



10" 



,3/2 



and 



logiCn®- 5 ) + Bn[ 

where the value of B « 0.47 cm 7/6 , C 
ni is 100 times the number density. Figure [2] illustrates the 
stiffening of the EOS as a function of number density. It can 
be seen from the figure that gas remains close to isother- 
mal up to 10 cm -3 and at higher densities opacity effects 
stiffen the value of 7 well beyond unity. An increase in the 
value of 7 enhances temperature well above 10 4 K. Even 
a value of 7 ~ 1 should already be sufficient to halt frag- 
mentation and will lead to isothermal collapse. Cooling due 
to helium becomes important at temperatures higher than 
5.0 xlO 4 K. Above 10 3 cm opacity effects due to thermal 
bremsstrahlung trapping come into play and cooling con- 
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Figure 2. Variation of gamma with density. The change in poly- 
tropic index (7 — 1) is plotted as a function of the number density 
(cm -3 ) of the gas. The figure shows the stiffening of the EOS by 
trapping Lyman alpha photons. 



tinues to be suppressed. A value of 7 ^ 4/3 leads to an 
adiabatic collapse. Helium opacity effects can be included 
b y making corrections to equation 2. Details can be found 
in ISpaans fe <j2Q06h . 



4 RESULTS 

Our simulations show that density perturbations decouple 
from the Hubble flow and start to collapse through gravi- 
tational instability. The fate of these perturbations depends 
on the cooling and dynamical time scales. At low densities 
cooling is very inefficient. We see in our simulations that ini- 
tially the gas is shock heated at low densities 10 -1 cm -3 ). 
Virialization processes transform gravitational potential en- 
ergy into kinetic energy of the gas and dark matter. As 
dark matter cannot dissipate, for gas this energy is con- 
verted into internal energy through shocks which raises the 
temperature of the gas to the virial temperature. We start 
with an adiabatic equation of state as perturbations in the 
inter-galactic medium are adiabatic in nature. According to 
linear perturbation theory, the EOS s hould show a power- 
law behavior in the low den sity regime (|Hui Gnedinlll997l : 
ISoringel & Hernauistll20 02). given by the equation 



T (p/po 



(3) 



where To and po are the cosmic mean density and temper- 
ature. Our results agree with this analytical temperature 
density relation for the intergalactic medium. 

At redshift 15, gas begins to cool and collapses in the 
gravitational potential provided by the dark matter. At den- 
sities of the order of 10 4 times the mean density of the Uni- 
verse cooling due to atomic lines (hydrogen and helium) 
becomes very efficient and gas becomes isothermal with a 
value of gamma of unity. Gas continues to collapse if the 
cooling time is smaller than the dynamical time. We have 
investigated the dynamic and thermal evolution of primor- 
dial gas which is cooled or heated and collapses under differ- 
ent post-shock heating mechanisms. In the next subsections 
the results are discussed for each case individually. We use 
the isovolume connected components finding algorithm im- 
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plem ented in the visualization program Visit ([Childs et al.l 
120051 ) to find fragments. For a given density threshold, we 
find the fragments and calculates their masses. 



4.1 Atomic line cooling 

Metal-free halos with virial temperatures > 10 4 K and no 
Lyman alpha trapping can be prone to molecular hydrogen 
formation if not irradiated by an intense UV flux. There 
may exist some halo pairs where molecular hydrogen is pho- 
todissociated by Lym an Werner radiation of t he other halo 



odissociated by .Lym an Werner radiation or t ne otner nalo 
Diikstra et al.l 12008") . Recent calculations by IShang et al.l 



( 2010h indicate that specific intensities required to suppress 
the formation of molecular hydrogen are 3-10 times lower 
than previous estimates of J cr i t ~ 1000 in units of 10 -21 erg 
cm -2 s -1 Hz -1 sr -1 . Consequently, the population of such ha- 
los is enhanced by a factor of « 10 3 . We assume that the halo 
is metal free and illuminated by intense UV flux. Hence, the 
formation of H2 does not take place. Initially gas heats up to 
its virial temperature. When the cooling time becomes less 
than dynamical time, it begins to cool and collapse in dark 
matter potentials. As cooling due to Lyman alpha photons 
is very efficient, the gas cannot virialize by gaining internal 
energy, so its kinetic energy must be increased to reach virial 
equilibrium. Consequently, during virialization the gas be- 
comes turbulent. It continues to cool and collapse into the 
center of a halo. 

Gas in the halo collapses up to densities of a few 
times 10 5 cm -3 . It fragments into two massive clumps of 
2.6 x 1O 6 M and 1.3 x 10 6 M o as shown in figure[3] Collapse 
is nearly isothermal and gas in the halo is cooled down to 
only a small fraction of the virial temperature. Isothermal 
behavior is depicted in temperature slice in figure [4] The 
temperature remains around ~ 8000 K. The value of the 
poly tropic index remains close to one. Isothermal collapse is 
preventing it to fragment into smaller clumps and the Jeans 
mass is 2 x 1O 5 M . 

The gas passes through different phases during collapse 
which are shown in the temperature-density phase diagram 
in figure [5] It is initially shock heated at low densities and 
subsequently cooled by atomic line cooling. Most of its mass 
is lying at low densities and only a small fraction (< 1%) is 
cooled and collapsed into clumps. It falls into the center of a 
halo through filaments and the direction of flow is shown by 
velocity vectors in figure [6] The typical velocity is around 
100 km We also examined the entropy of gas in the 
halo and found that it is 10 -2 keVcm -2 . 

Our results are in agreement wi th previous results 
(jBromm fe Loebll2003l : IWise et alJl2008h . We do not see disk 
formati on in our case as mentioned by iRegan fc Haehneltl 
(2009b). The reason for this disagreement could be due to 
difference in resolution. We do not exclude the possibility 
of disk formation which may form after collapse. Inefficient 
fragmentation of a halo for isothermal collapse is according 
to the expectation of theoretical models. 



4.2 Molecular hydrogen cooling 

Adding molecular hydrogen cooling enhances the low tem- 
perature cooling and increases the fragmentation. We see 
that gas is heated up to its virial temperature by shock 



waves at lower densities and subsequently cools by atomic 
line cooling. At densities higher than 1 cm -3 cooling due to 
molecular hydrogen comes into play and lowers the temper- 
ature down to rsj 100 K. During collapse it passes through 
different phases, its multiphase behavior is shown in the 
temperature-density phase diagram in figure Most of the 
gas is lying at low densities and only a small fraction (< 1%) 
is cooled and collapses into minihalos. 

The equation of state is softened due to H2 cooling and 
the value of the polytropic index is lowered to less than unity. 
Hence, the Jeans mass is reduced. A density slice is shown in 
figure[8] H2 cooled gas is seen at the intersection of filaments. 
The gas collapses into 6 minihalos in contrast to the Lyman 
alpha trapping case, where gas does not collapse into mini- 
halos. Their masses are 2 x 1O 6 M , 8 x 1O 5 M , 8 x 1O 5 M , 3 x 
1O 5 M , 1 x 1O 5 M and 3.4 x 1O 5 M . Only two of the mini- 
halos are visible in the slice. Their sizes are significantly 
smaller than the other cooling cases. The Jeans mass is re- 
duced to <~ 1O 3 M and the gas collapses to densities of a few 
times 10 6 cm -3 . The gas becomes turbulent during virializa- 
tion and f alls into the center of a halo through cold streams 
as seen in (IKeres et al. 2009; Fardal et al. 2001; Dekel et al.l 



120091 : iGreif et alJ l2QQ8: Wise & Abel 2003). The density is 



10 cm - and temperature is <~ 300K in these streams. 
The direction of flow is presented in figure [9] The tempera- 
ture of the minihalos is around 100 K and is shown in figure 
1101 Shock heated gas is seen in the vicinity of the minihalos 
and the penetration of cold flows is even more prominent in 
this plot. The introduction of H2 cooling has significantly 
reduced the entropy of the gas down to 10 -5 keVcm -2 . En- 
tropy is low in the minihalos, higher in outskirts, and even 
boosted in shocks (10 2 keVcm -2 ). 



4.3 Lyman alpha trapping 

It is generally assumed that gas, even though very optically 
thick (r > 10 7 ) to scattering of Lyman alpha photons, can 
cool efficiently through Lyman alpha because radiation leaks 
out in the line wings. High columns of neutral atomic gas 
make the gas so optically thick that photon travel time ex- 
ceeds the tff. So, using the effectively optically thin cooling 
rate is not a good approximation above 10 20 cm -2 . Conse- 
quently, the EOS stiffens as explained in the previous sec- 
tion and cooling is suppressed. Gas remains almost optically 
thin up to densities on the order of 10 cm -3 . At densities 
^ 10 2 cm -3 it becomes opaque enough that cooling due 
to Lyman alpha photons is quenched. When the gas den- 
sity exceeds 10 3 cm -3 , at typical columns of the order of 
10 23 cm -2 , the 2-photon continuum channel is shut down 
and photo ns are absorbed by n eutral hydrogen in the sur- 
roundings ([Spaans fc Silkll2006h . Under these conditions the 
halo collapses adiabatically to form a massive object. 

Gas is initially shock heated at low densities, becomes 
isothermal at densities ~ 1cm -3 , while at densities ^ 1cm -3 
EOS is stiffened because of Lyman alpha trapping. This 
multi-phase behavior is depicted in the density-temperature 
phase diagram in figure ITT1 Most of the gas is lying at low 
densities and only a small fraction (< 1%) collapses to high 
densities. Gas in the halo becomes turbulent during virial- 
ization and is segregated according to its angular momentum 
distribution (low angular momentum gas falls into the cen- 
ter). Hence, the collapse happens in the center of a halo. An 
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Figure 3. Density slice through through the center of a halo 
for the atomic cooling case. Values corresponding to each color 
are shown in color bar, where blue color represents the low- 
est density and red represents the highest density. Distance 
scales are in comoving units. Here lkpc in comoving units 
corresponds to 106 pc in proper units. 
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Figure 5. Density-temperature phase diagram for the atomic 
cooling case. Temperature is plotted against volume weighted 
density in filled color contours by merging adaptive grid into 
uniform grid, where red color shows the highest mass and pur- 
ple shows lowest mass. 
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Figure 4. Temperature slice through the center of a halo for 
the atomic cooling case. The figure shows the temperature cor- 
responding to density slice of figure [3] Blue color represents 
the lowest temperature and red represents highest tempera- 
ture. 




Figure 6. Velocity vectors are overplotted on density slice 
for the atomic cooling case. Velocity vectors are indicated by 
white color arrows. Head of the arrows indicates the direction 
of flow. 



object of mass 4.0 x 1O 6 M0 is formed at the intersection of 
filaments which may lead to a massive object i.e. an IMBH. 
We see signs of a cold flow, which has already been seen 
in numerical simulation s by (|Dekel et al.ll2009l ; iKeres et al.l 
I2009I : iFardal et al.ll200lh . Baryons accumulate into the cen- 
ter of a halo via a two-phase medium by entrance of cold 
flow streams through the shock heated medium. These flows 
also explain the observ ed Lyman alpha blob morphologies 
(iDiikstra fe Loebll2009h . The direction of flow is shown by 
velocity vectors in figure 1 121 The values of the velocity in 
the center of halo are about 100 km s _1 which are typical 
values for such objects. 

The morphology of the object is shown in figure IT3l Dif- 
fuse hot gas is visible in the vicinity of the object. Contrac- 
tion of gas is halted and the Jeans mass is boosted due to Ly- 
man alpha trapping. Diffuse gas may feed the central object 



during later stages of collapse. The halo does not collapse to 
high densities (few xl0 3 cm -3 in this case) as compared to 
other cases because of suppressed cooling. The temperature 
is enhanced due to trapping of Lyman alpha photons as de- 
picted in figure [141 It is ^ 10 4 K in the center of the halo and 
even higher in the shock fronts. Trapping of Lyman alpha 
photons keeps the gas H2 free, as collisional dissociation and 
charge exc hange at a temperatu r e above 5000 K destro ys H2 
formation (jOh fc Haimanll2002h . fSoaans & Silkl (|2006l ) have 
confirmed through models that column densities in the range 
of 10 22 - 10 25 cm will be sufficient to self shield the gas 
even from X-rays. This puts our arguments in favor of sup- 
pressed fragmentation on sounder footing. The inability of 
gas to fragment makes such halos prone to the formation of 
massive objects. The collapse is never halted but only de- 
layed in the case of Lyman alpha trapping as compared to 
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Figure 7. Density-temperature phase diagram for the molec- 
ular cooling case. Temperature is plotted against the volume 
weighted density in filled color contours by merging adaptive 
grid into a uniform grid, where red color shows the highest 
mass and purple shows lowest mass. 





Figure 8. Density slice through the center of a halo for the 
molecular cooling case. Values corresponding to each color are 
shown in the color bar, where blue color represents the low- 
est density and red represents the highest density. Distance 
scales are in comoving units. Here lkpc in comoving units 
corresponds to 106 pc in proper units. 




Figure 9. Velocity vectors are overplotted on the density 
slice for the molecular hydrogen cooling case. Velocity vectors 
are indicated by white color arrows. The head of the arrows 
indicates the direction of flow. 



Figure 10. Temperature slice through the center of a halo 
for the molecular hydrogen cooling case. The figure shows the 
temperature corresponding to the density slice of figure [8] The 
temperature of fragments is close to 100 K which can be seen 
in blue color in the plot. 



other simulations (molecular hydrogen and atomic cooling). 
This delay is caused by the enhanced radiation pressure as 
more work needs to be done on the gas. Entropy of gas is 
increased as a consequence of Lyman alpha trapping up to 
0.1 keVcm 2 . It is even higher (10 keVcm 2 ) in shocks. We 
see that trapping of Lyman alpha photons increases the en- 
tropy of gas which is close to the lower limit for entropy 
in simul ations of galaxy grou ps and clusters with radiative 
cooling ([Borgani et al.ll2005h . We do not include feedback 
effects which may inject extra energy and raise the entropy 
of the gas to even higher values. 



5 DISCUSSION AND CONCLUSIONS 

We have probed the influence of thermodynamics on the 
formation of the first objects in primordial gas. We have in- 
troduced the impact of Lyman alpha trapping, atomic and 
molecular hydrogen cooling, and performed cosmological 
three-dimensional simulations down to a redshift of 8.4 using 
the AMR code FLASH. We exploited the AMR technique 
to add dynamic levels of refinement wherever required, so 
that our resolution is sufficient to resolve the most essential 
physics. We have assumed that halos are metal- free. Pockets 
of m etal-free gas can exist until redshift ~ 6 (| Trent i et al.l 
l2009h . 

Line trapping of Lyman alpha photons stiffens the EOS 
and enhances the Jeans mass. We noticed that Lyman al- 
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Figure 11. Density-temperature phase diagram for the Ly- 
man alpha trapping case. Temperature is plotted against vol- 
ume weighted density in filled color contours by merging adap- 
tive grid into uniform grid, where red color shows the highest 
mass and purple shows the lowest mass. 
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Figure 13. Density slice through the center of a halo for the 
Lyman alpha trapping case. Density is shown in rainbow colors 
where red represents highest density and purple represents 
lowest density. Distance scales are in comoving units. Here 
lkpc in comoving units corresponds to 106 pc in proper units. 
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Figure 12. Velocity vectors are overplotted on density slice 
for the Lyman alpha trapping case. Velocity vectors are indi- 
cated by white color arrows. Head of the arrows indicates the 
direction of flow. 
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Figure 14. Temperature slice through the center of a halo for 
the Lyman alpha trapping case. The figure shows the temper- 
ature corresponding to density slice of figure [13] Temperature 
of object is high compared to atomic cooling case. Distance 
scales are in comoving units. 



pha photons trapped in halos inhibit cooling below ~ 10 4 K, 
consequently, gas fragmentation is halted and a massive ob- 
ject is formed at the center of a halo. We performed a set 
of simulation with one higher level of refinement for this 
case and find similar results as mentioned before. As trap- 
ping continues, formation of H2 will be suppressed and a 
halo will collapse adiabatically. The outcome of the object 
depends on how fast and efficient it accretes mass. It can 
directly collapse into a black hole or a supermassive star 
as an interim stage. We conclude that trapping of Lyman 
alpha photons during gas cloud collapse provides a venue 
for black holes formation, particularly when there is modest 
background UV flux to suppress molecular hydrogen cool- 
ing. Radiative transfer effects thus play a crucial role in the 
formation of massive objects observed at z > 6. 

Our results are only valid for metal free halos. The 
addition of trace amounts of metals can change the sce- 



nario for Lyman alpha trapping. If dust is present then the 
gas temperature is lowered, and the equation of state soft- 
ened, due to absorption and IR re-emission of Lyman alpha 
photons. Halos might have been pre-enriched by the heavy 
metals formed in and d i spersed through sup ernovae explo- 
sions ([Greif et al.l l201ol : Ijohnson et al. [2008). Such a sce- 
nario would suppress the effectiveness of Lyman alpha trap- 
ping because dust would efficiently absorb the Lyman alpha 
photons and re-radiate them in the far-infrared. As such, we 
again find an upper limit on the importance of Lyman al- 
pha trapping. Still, if the medium is inhomogeneous, which 
it very likely is, then Lyman alpha may scatter of dense gas 
condensations without being absorbed b y dust and thus re- 
main trapped in a collapsing structure (jHaiman &; Spaansl 
Il999l : iNeu feld 1991). Fragmentatio n of the gas is inevi table 
in the presen ce of dust and metals (|Omukai et al.l l2008). Re- 
cent work bv lSchleicher et all (|201(F shows that after initial 
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Figure 15. The figure represents a higher density peak halo for the Lyman alpha trapping case. The upper left panel shows the density 
slice through the center of a halo. The temperature slice corresponding to density slice is depicted in the upper right panel. The bottom 
left panel shows the phase diagram. The flow of gas is seen in the bottom right panel. Values corresponding to colors are shown in the 
color bars. Distance scales are in comoving units. Here lkpc in comoving units corresponds to 53 pc in proper units. 



collapse, cooling may proceed through other atomic states 
in the presence of Lyman alpha trapping. It can soften the 
EOS to unity and lead to disk formation, which may or may 
not fragment to form a top-heavy stellar initial mass func- 
tion. We have explored the important phase of Lyman al- 
pha trapped halos. Cosmological simulations should be per- 
formed to examine the effect of other cooling channels like 
2s- Is and 3-2 transitions in halos that are optically thick to 
Lyman alpha photons. 

For the atomic cooling case, we find that the halo col- 
lapses into one or two massive clumps. As long as H2 cooling 
is suppressed, gas collapses isothermally and remains hot 
having temperature around 8000 K. Further fragmentation 
is unlikely. Therefore, such halos are plausible candidates 
for the formation of intermediate mass black holes or su- 
permassive stars as an intermediate stage. The formation of 
molecular hydrogen does not take place if internal or exter- 
nal UV radiati on from the star forming halo can sup press 
H 2 formation (jShang et alJl20ld : [Piikstra et alJ l2Q08h . At 
densities > 10 8 c m~ 3 molecular hy drogen may form due to 
3-body reactions ( Pall a et al .11 19831 ) in the absence of intense 
UV flux and may induce fragmentation. The further fate of 



these objects should be explored through cosmological sim- 
ulations by including detailed physics. 

We have also performed simulations by selecting a high 
density peak halo. We noticed that in the higher density 
peak halo the free-fall time becomes shorter and a halo col- 
lapses about redshift 25. We followed the collapse down to 
redshift 18 and compared the results of atomic line cooling 
with Lyman alpha trapping case. We find that the Lyman 
alpha trapped halo does not fragment irrespective of the 
high or low density peak in contrast to the atomic cooling 
case. Trapping of Lyman alpha photons keeps the gas hot 
(T^ 10 4 K) and contraction of a halo is halted. The Jeans 
mass becomes higher and it does not fragment. The size of 
the fragment is larger than the atomic cooling case. The 
mass in a higher density peak halo is 5 x 1O 7 M0. The mass 
of the object is not much higher than the atomic cooling case 
because more gas falls into the halo due to efficient atomic 
cooling and condenses. The density, temperature, phase di- 
agram and velocity vectors for a higher density peak halo 
are shown in figure ITol 

For the atomic cooling case, morphology and thermody- 
namical properties of the clumps are similar in both higher 
and lower density peak halos. The density, temperature, 
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Figure 16. The figure represents a higher density peak halo for the atomic cooling case. The upper left panel shows the density slice 
through the center of a halo. The temperature slice corresponding to the density slice is depicted in the upper right panel. The bottom 
left panel shows the phase diagram. The flow of gas is seen in bottom right panel. Values corresponding to colors are shown in the color 
bars. Distance scales are in comoving units. Here lkpc in comoving units corresponds to 53 pc in proper units. 




Figure 17. Density slices of the inner 2 kpc region (zoomed) of the figures \W\ (upper left panel) and figure [15] (upper left panel). The 
left panel shows the density slice for the atomic cooling case. The right panel shows the density slice for the Lyman alpha case. Contours 
are overplotted to show the density structure. Distance scales are in comoving units. Here lkpc in comoving units corresponds to 53 pc 
in proper units. 
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Figure 18. The figure represents a higher density peak halo for the molecular hydrogen cooling case. The upper left panel shows the 
density slice through the center of a halo. The temperature slice corresponding to the density slice is depicted in the upper right panel. 
The bottom left panel shows the phase diagram. The flow of gas is seen in the bottom right panel. Values corresponding to colors are 
shown in color the bars. Distance scales are in comoving units. Here lkpc in comoving units corresponds to 53 pc in proper units. 




Figure 19. Velocity field corresponding to [171 The left panel shows the velocity field for the atomic cooling case. The right panel shows 
velocity field for the Lyman alpha case. Velocity vectors represent the flow of gas after subtracting back ground flow. 



12 M.A. Latif et al. 



phase diagram and velocity vectors for a higher density peak 
halo are shown in figure 1161 We find that the higher den- 
sity peak halo fragments as in lower density peak case. It is 
depicted in density slice (zoomed in) in figure rm Density 
contours are overplotted on density slice and are compared 
with Lyman alpha case. One can see that there is more struc- 
ture in atomic case as compared to trapping case. We found 
that clump surrounding massive one is gravitationally bound 
and mass is close to the Jeans mass. We also noticed that 
smaller clump is rotating around massive one. The veloc- 
ity field for the atomic cooling case (inner 2 kpc region) is 
shown in figure [19] and is compared with Lyman alpha case. 
There is also divergence in the flow of gas corresponding 
to density contours. It supports our claim of fragmentation. 
The masses of the clumps in higher density peak halo are 
- 3 x 1O 7 M ,4.O x 1O 5 M . These fragments will become 
even more prominent in higher resolution simulations. 

In the case of molecular hydrogen cooling, the gas is 
cooled down to 100 K and collapses into minihalos. The 
minihalos formed via H2 cooling are the potential sites for 
Pop III star formation. Following the collapse becomes a 
computationally difficult task as the Jeans mass keeps drop- 
ping. We stop the simulation when the Truelove criterion 
is about to be violated. Adding higher resolution may help 
in fragmentation. Further fragmentation may occur at the 
higher densities due to an enhanced fraction of the molecular 
hydrogen formed through 3-body reactions. We did not in- 
clude the HD cooling which is a very efficient coolant around 
100K. Introduction of the HD cooling will cool the gas below 
100K and may enhance fragmentation. 

We have performed a set of simulations for the higher 
density peak H2 cooled halo and compared it with the lower 
density peak halo. We find that the gas collapses into mini- 
halos. Efficient fragmentation of H2 cooled gas halo occurs 
regardless of being in low or high density peak halos. The 
thermodynamical properties of the minihalos are very simi- 
lar in both cases. The difference in the collapse time of the 
higher and lower density peak halos is the same as mentioned 
for the atomic cooling case. The density, temperature, phase 
diagram and velocity vectors for a higher density peak halo 
are shown in figure [151 The masses of the minihalos in this 
case are 9 x 1O 5 M , 6 x 1O 5 M , 5 x 1O 5 M , 3 x 1O 5 M , 2.6 x 
1O 4 M , and 9 x 1O 4 M . Stars are most likely to form in the 
H2 cooled minihalos. The sizes of the clumps are bigger than 
H2 cooled fragments due to difference in the Jeans masses. 
Fragmentation to lower than the Jeans scale is not possible 
as gas becomes pressure supported. 

We also see that the entropy of the gas is increased up 
to 0.1 keVcm 2 by trapping Lyman alpha photons. It is in- 
creased in shocks and lost due to cooling. Complete entropy 
calculations are beyond the scope of this article. We have 
not included the feedback effects from the first objects in 
our simulations. The presence of UV flux from Pop HI stars 
will suppress the formation of H2 (|Diikstra fe Loebl 120091 : 
IShang et alJl2010h . Further fragmentation will be inhibited 
in these halos. It will support the mode of massive objects 
formation in the halos that we consider here. It should be 
noted in this that star formation is not completely sup- 
pressed by a UV background due to efficien t self-shielding 
or positive feedback that results from shocks ([ Johnson et al I 
120071 : iRicotti et al. l l2002h . X-rays from miniquasars can stim- 
ulate the formation of molecular hydrogen (|Haiman et al.l 



l2000h . High columns present in our case will partly self shield 
the gas from X-rays feedback but the gas in the surround- 
ing might be effected. Removal or redistribution of angular 
momentum is an obstacle in the formation of massive ob- 
jects. The gas can transfer angular momentum via shoc k 
waves and bar instabilities (IWise et alJl2008l : lLarsonll2009h . 
We noticed that gas becomes turbulent during virialization 
and helps in the dissipation of angular momentum. We do 
see strong shock waves occurring at low densities in our sim- 
ulations which also play important role in the distribution of 
angular momentum. Numerical simulations should be per- 
formed to get detailed insight about angular momentum is- 
sue. 

In any case, thermodynamical properties of the gas play 
an indispensable role in the fragmentation of gas cloud and 
the formation of primordial objects. Our results include ra- 
diative transfer effects of the Lyman alpha photons in an 
approximate manner, but still give good insight into radia- 
tive transfer effects on primordial gas physics. Forthcoming 
telescopes like LOFAR and JWST will be able to observe 
the objects formed in the early Universe and will help to 
scrutinize the present paradigm of primordial objects for- 
mation. 
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